[T E E2 O] = textread('../data/spin.dat', '%f,%f,%f,%f');
plot(T,E,'r');
axis([0 300 -300 300]);
